Light bosons and photospheric solutions to the solar abundance problem 
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It is well known that current spectroscopic determinations of the chemical composition of the 
Sun are starkly at odds with the metallicity implied by helioseismology. We investigate whether 
the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. 
We examine the impact of particles with axion-like interactions with the photon on the inferred 
photospheric abundances, showing that resonant axion-photon conversion is not possible in the 
region of the solar atmosphere in which line-formation occurs. Although non-resonant conversion in 
the line-forming regions can in principle impact derived abundances, constraints from axion-photon 
conversion experiments rule out the couplings necessary for these effects to be detectable. We show 
that this extends to hidden photons and chameleons (which would exhibit similar phenomenological 
behaviour), ruling out known theories of new light bosons as photospheric solutions to the solar 
abundance problem. 



I. INTRODUCTION 

Agreement between predictions of solar interior models 
and helioscismological measurements has deteriorated in 
recent years, due to a downwards revision in standard so- 
lar abundances. Recent photospheric analyses [1-9] have 
indicated that the solar metallicity is about 20% lower 
than previously thought [10]. In particular, the current 
reference photospheric abundances [9] of the important 
elements C, N and O are smaller than the previous refer- 
ence values [10] by 0.09, 0.09 and 0.14 dex, respectively. 
Standard solar models computed with the revised sur- 
face metal abundances show large discrepancies with the 
depth of the convection zone, the surface helium abun- 
dance and the sound speed profile inferred from helio- 
seismology [11-16]. This has been variously referred to as 
'the solar model(ling) problem' [17-19], 'the solar oxygen 
crisis' [20-22], or 'the solar abundance problem' [23, 24]. 

The most common approaches to the solar abundance 
problem to date have been to recheck the photospheric 
abundances [18-21, 25-27], or to modify solar models in 
an attempt to accommodate both the measured abun- 
dances and data from helioseismology. Efforts in the lat- 
ter category have included enhanced diffusion [28], late- 
stage accretion of low-mctallicity gas [23, 24, 28, 29], 
modifications of opacities [30-32] , internal gravity waves 
[33, 34] and enhanced energy transport in the solar core 
due to gravitational capture and scattering of dark mat- 
ter [35-37]. Despite nearly a decade of effort (and a few 
errant claims to the contrary), none of these approaches 
has proven successful in solving the problem. 
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Here we propose an alternative approach: examining 
effects of new physics on the measurement of solar photo- 
spheric abundances. Specifically, we investigate the im- 
pact on derived elemental abundances of mixing between 
photons and new light bosons, such as axions, in the line- 
forming regions of the photosphere. 1 

The 5770 K thermal spectrum of the Sun is formed 
around optical depth r = 1, which defines the location of 
the solar 'surface'. For a height z above the surface, 

/>oo 

tx(z) = / a x (z')dz'. (1) 

J z 

Here ct\(z) = K c x (z)p(z) is the linear attenuation coeffi- 
cient at wavelength A, where k c x is the continuum opacity 
and p is the local gas density. Solar abundances are in- 
ferred from absorption lines. These are formed as atomic 
or molecular bound-bound transitions increase the opac- 
ity in small wavelength ranges around the transitions. 
The higher opacity in these spectral lines means that the 
light at those wavelengths escapes from layers above the 
continuum-forming layers (the photosphere). As the gas 
temperature decreases with height, the line-centre pho- 
tons escape from cooler layers than the continuum pho- 
tons, giving rise to absorption lines. 

The equivalent width of a line (the integrated area be- 
tween the line and continuum) is a relative measure of 
line strength, and depends not only on the abundance of 
a given element, but also on the attenuation of the total 



A more complicated but related scenario was considered in [38], 
where high-energy axions produced in the solar core would recon- 
vert to photons in the photosphere and exert radiation pressure 
preferentially on heavier elements, driving them from the photo- 
sphere. 
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continuum as it passes through the line-forming region. 
For this reason, both line and continuum fluxes must be 
carefully calculated with radiative transfer simulations 
before elemental abundances can be inferred. If addi- 
tional plasma effects in the photosphere enhance contin- 
uum photon attenuation via oscillation into light bosons 
such as axions, then the equivalent width of absorption 
lines used to calculate elemental abundances would be 
suppressed. By combining the axion conversion effect 
with earlier radiative transfer calculations, we will show 
how this could occur in principle, and assess whether this 
effect could arise from existing particle theories in such 
a way as to solve the solar abundance problem. 

Historically, the axion field arose as a candidate solu- 
tion to the strong CP problem of QCD. Although the 
term 

e 3 G^G^, (2) 

is not forbidden in any way in the QCD Lagrangian, it 
seems to be absent from the standard model (SM) of 
particle physics. Here G^ v is the gluon field strength, 
and its dual is G^ u = e a ^^ u G a ^. The coupling constant 
03 is experimentally consistent with zero, so some sort 
of dynamical method is required in order to force it to 
be small without fine-tuning. The Peccei-Quinn (PQ; 
[39]) mechanism introduces a broken U(l) symmetry that 
naturally yields this cancellation, with the added effect of 
producing a new light, propagating, pseudoscalar degree 
of freedom: the axion. 

The proposed existence of the axion field has given rise 
to a rich variety of phenomenological predictions. These 
include axion production in the dense cores of stars ([40] 
and references therein), "light shining through wall" ef- 
fects [41], and many theories of axionic dark matter ([42] 
and references therein). Many of these effects are due 
to the effective coupling £ a77 that must arise between 
the PQ axion field a and the photon gauge field strength 

C ail = ^aF^F^ = -g ai aE ■ B (3) 

with coupling constant g ai . It is this axion-photon- 
photon vertex that allows axions to convert into photons 
in the presence of electric (E) or magnetic (B) fields, 
and vice versa. The mass and coupling of a PQ axion 
must be very small (with m a oc l/g ai ), and currently 
remain (mostly) far below the threshold of experimental 
searches. A rough lower limit of 10 _6 eV can but put 
on the axion mass, in order to avoid overdosing the uni- 
verse, and an upper limit of m a < 10 -3 eV ensures that 
stellar evolution is not significantly modified [43]. 

The PQ mechanism in QCD is not the only possible 
origin of such a particle, however. Axion-like particles 
(ALPs) are a ubiquitous element in string theory com- 
pactifications [44, 45], in which the moduli describing the 
compact extra dimensions can appear to observers as in- 
ternal U(\) gauge degrees of freedom. More generally, 
U(l) symmetries appear in many high-energy extensions 



of the SM, giving rise to exotic particles with interactions 
similar to the axion. Like axions, ALPs interact with 
photons via a term of the form Eq. 3. Unlike axions, 
they do not necessarily solve the strong CP problem, nor 
require such small couplings or m a oc l/g ai . 

We begin by showing in Section II how photon- ALP 
oscillation can affect the linear attenuation coefficient - 
and hence the equivalent width of an absorption line - 
in the solar atmosphere at optical depths r < 1. We 
demonstrate in Section III that with a large enough cou- 
pling g ai , a low-mass particle with a term like Eq. 3 could 
affect the solar atmosphere enough to modify absorption 
lines in an observable way. As we discuss in Section IV 
however, the coupling strength required to produce such 
an effect is much larger than the coupling allowed by cur- 
rent experimental bounds on standard ALPs. We go on 
to show that this rules out not only ALPs as photospheric 
solutions to the solar abundance problem, but also hid- 
den photons and chameleons. We conclude in Section V 
that unless some other light boson with an axion-like or 
similar kinetic-mixing type coupling to the photon can 
be found that somehow evades the existing experimen- 
tal constraints, the solar abundance problem cannot be 
solved by light bosons in the photosphere. 

II. THEORY 

A. Changing inferred abundances 

We first consider the standard case, without ALPs. 
On the linear portion of the curve of growth, the equiv- 
alent width of an absorption line W\ is approximately 
proportional to the local ratio of the line-centre (k\) to 
continuum (k^) opacities at the height of line formation 
[46]: 

W x cx^= Vx . (4) 

The abundance-dependence of W\ is encoded in -q\ by 
way of k\, which is itself proportional to the number of 
absorbers along the line of sight. The number of ab- 
sorbers is given by the number of atoms with electrons 
in the lower level of the transition, which is some known 
(temperature- and density-dependent) fraction of the to- 
tal abundance of the element in question, e. 

In the presence of ALPs, the local continuum opac- 
ity receives an additional effective contribution from 
photon- ALP conversion. In this case 

W\ oc — ^ = -9*-, (5) 

where we have defined the 'axionic extinction ratio' rf^ = 
k a/ k a m ana l°gy with the line extinction ratio r]\. The 
abundance-dependence of Eq. 5 is identical to that of 
Eq. 4, i.e. linear in e, with the same constant of propor- 
tionality. However, given that W\ is known from obser- 
vation, analysis using either Eq. 4 or Eq. 5 is constrained 
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to return the same answer for W\. Denoting the abun- 
dance derived by ignoring ALPs (Eq. 4) as e st d, and that 
including the effects of ALPs (Eq. 5) as e COIT , we have 



W x (Eq. 4) = 6 std 
W A (Eq. 5) e corr /(l + Tfi) 



= 1. 



(6) 



The abundance correction to be applied to the standard 
result in order to obtain the ALP-corrected abundance is 
thus 

A log e = log e corr - log e std = log(l + r}*). (7) where: 

With Eq. 7 in hand, the remaining task is then to calcu- 
late rf^ in the line-forming regions of the photosphere. 



B. Photon-ALP conversion in the Solar atmosphere 

Regardless of the precise ALP model, the presence of 
a coupling term with the photon (Eq. 3) implies that 
under the right conditions, a thermal photon in the pho- 
tosphere may oscillate into an ALP and vice-versa. This 
occurs in the presence of an external photon field: either 
a bulk magnetic field, or the oscillating fields of the stellar 
plasma. Only the former effect is relevant in the region of 
interest for line formation, as plasma field-induced tran- 
sition rates (Eq. 9) turn out to be extremely small in the 
photosphere. In order to arrive at the quantity 77° we first 
compute the linear attenuation coefficient a^, related to 
the ALP opacity by: 



(8) 



Conversion in plasma fields 



The transition from photon 7 to ALP a induced by 
the oscillating electric fields in a plasma occurs mainly 
from resonant Primakoff conversion. Neglecting recoils, 
the transition rate is given by [47, 48] : 
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where k s is the screening scale in the Debye-Hiickel ap- 
proximation: 



(10) 



Zj is the charge of the i th ion and the number densities 
n c , n; of electrons and ions, respectively, are found by 
solving the equation of state. This effect is highly sup- 
pressed at the relatively low charge densities of the solar 
atmosphere, and contributes a negligible amount at op- 
tical depths of interest to us. 



2. Conversion in a bulk magnetic field 

Conversion in a bulk magnetic field can be much more 
significant than in the plasma field. The probability that 
a photon polarised in the direction of the B field will 
have converted to an ALP after travelling a distance Sz = 
z — zq from height zq to height z is [49, 50]: 
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Bt is the component of the magnetic field transverse to 
the direction in which the photon propagates, u) = 2itc/\ 
is the frequency of the incoming photon, and the plasma 
frequency is u> p \ — 4nan c /m c , where a is the fine- 
structure or electromagnetic coupling constant. Due to 
the E ■ B form of the interaction (3), only the photons 
with polarisation parallel to the magnetic field are sus- 
ceptible to conversion; this is the meaning of the 1 1 sub- 
script. Ajp = w(n[] ac - 1), where n[] ac ~ 1 + {aB/4nB c ) 
is the "vacuum" contribution to the refractive index in- 
duced by the magnetic field. The critical field B c — 
m 2 /e ~ 10 12 eV 2 ~ 4x 10 9 T is much larger than the 
fields B ~ 0.1 T that we consider in this work; A|j ac is 
therefore negligible for all cases of interest here. Simi- 
larly, A^ as = w(n^ as — 1), where n^ as is the contribution 
to the refractive index by neutral hydrogen gas. A rela- 
tivistic calculation of the static polarizability of hydrogen 
[51] gives a (Hi) = 6.6679437 x 10~ 25 cm 3 . The index of 
refraction is then: 



,gas 
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= a/1 +47rn(Hl)a (Hl) 



(13) 



where n(Hi) is the number density of neutral hydrogen 
atoms. This is valid in the low densities of the solar 
atmosphere and as long as we are sufficiently far from 
the Lyman a line. 

Eq. 11 describes a Lorentzian (Breit-Wigner) distri- 
bution in A a oc 1/uj modulated by the sin 2 -term, with 
resonance at the plasma frequency and width A a7 cx B. 
When A|| — A a goes to zero, photon conversion to ALPs is 
enhanced. We later show that for m a ~ 10~ 6 -10 -3 eV, 
this has the potential to occur in the region of the so- 
lar atmosphere where line- format ion occurs. The plasma 
frequency uo v \ increases quickly with increasing optical 
depth. Given the w~j 4 suppression of Eq. 11, this means 
that conversion is highly suppressed at depths below the 
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FIG. 1. Green: average plasma frequency u p i (eV) as a 
function of optical depth at 500 nm (bottom axis) and av- 
eraged height (upper axis) in a convection simulation of the 
Solar atmosphere. Red: contribution of the neutral hydrogen 
(wAy as )'' 2 at photon energies corresponding to A = 300 nm 
(dashed) and 900 nm (solid). When the difference between 
the plasma and neutral hydrogen quantities (thick blue lines) 
is positive, oJ^y — cjA^ as > 0, resonant enhancement may occur. 
The value of this difference at a given height corresponds to 
the ALP mass m a for which resonant conversion takes place. 



resonance depth h les . There are thus three qualitatively 
different regions in the solar atmosphere: 1) Below h res , 
conversion is suppressed by k>~, 4 , ensuring that very little 
conversion occurs at depth; 2) at /i ros , resonance occurs: 
conversion is enhanced and Pq —¥ 1 for a small region; 
and 3) above h Tes , Pq becomes approximately constant, 
providing moderate conversion rates. The existence of 
region 2) is contingent on the plasma effects being larger 
than the contribution from neutral gas; otherwise, A|| 
and — A a are always positive, prohibiting cancellation. 

We point out that many references (e.g. [47, 50]) do not 
include A^ as . This is because in situations such as hot 
stellar interiors or the interstellar medium, the ionisation 
fraction is so large that A^ as is always subdominant to 
A p i. Indeed, the intermediate temperatures encountered 
in the solar atmosphere make it one of the few places 
where both the gas refractive index and plasma frequency 
contribute to similar degrees. Fig. I illustrates these ef- 
fects, showing that resonant conversion may only occur 
below logr ~ (i.e. at larger optical depths). The blue 
lines show the ALP mass at which resonant conversion 
may occur, for two different photon wavelengths. 

We now compute the effect of photon- ALP conversion 
on the photon intensity, for two complementary magnetic 
field configurations, starting with the case where the size 
of the magnetic domains s is larger than the scale height 
over which line formation occurs. After propagation over 
a distance Sz, photon loss to photon- axion conversion 



gives the change in intensity [50]: 

I 



I 7 (z) 



1 - -P Q (Sz) ) 7 7 (z ). 



(14) 



This assumes that the change of the photon intensity due 
to radiative transfer is small over the distance Sz and does 
not feed back on the calculation of the conversion rate. 
The linear attenuation coefficient due to ALP conversion 
is then 



Uz) 



A osc sin(A osc fe) 
- 1 + cos(A osc Sz) 



(15) 



where the prime symbol denotes differentiation with re- 
spect to z. When computing the full radiation transfer, 
(15) must then be added to the linear attenuation due to 
opacity. 

For a line with some mean optical depth of formation 
r, the range of heights over which will contribute to 
the line profile is approximately given by the scale height 
H K with which the continuum opacity changes by a factor 
of e, 
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This simply states the fact that to a rough approxima- 
tion, the width of a line's contribution function is given 
by the opacity scale height. 2 The relevant contribution 
of to the continuum opacity for any line formed at 
depth r can then be estimated by averaging over H K 
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where x(8z) = cos(A osc 5z). Integrating gives 



[1 - cos(A osc # K )] 



(18) 

Away from resonance, A osc 3> A a7 and the log may be 
Taylor-expanded about 1, giving 



1 



[1 - cos(A osc if K ) 



(19) 



Because we are interested in the resonance region as well, 
we will only make use of the full expression (Eq. 18). 

If photons instead propagate through n domains of 
length s with equal magnetic field strength but random 



2 More correctly, the scale height of the line-to-continuous opac- 
ity ratio rjx- As this scale height varies on a line-by-line basis, 
it is not especially useful for estimating effects on abundances 
as a function of formation height, in the transition-independent 
manner we do here. 
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orientation (such that Sz 
to [50]: 



ns 



) , their intensity is reduced 



(20) 



Here I a refers to the intensity of the ALP field, and Pq = 
Pq(s). If the ALP population is initially zero - we return 
to the validity of this assumption below - the photon 
intensity becomes : 
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If the opacity scale height H K (the propagation length of 
interest) is much larger than the magnetic domain size s, 
we can use the large n limit of Eq. 21 [50] 
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(22) 
(23) 



This is valid for a photon passing through a large num- 
ber Sz/s of magnetic domains of random orientation, but 
with identical field amplitudes, B. This serves the same 
function as the averaging procedure (Eq. 17). Note that 
Po is again evaluated at s, the size of the magnetic do- 
mains. The linear attenuation coefficient is then 



W = Po 3 

I 7 (z) 2,sl + 2exp(^) 



(24) 



Away from the resonance region, terms that are higher- 
order in Pq may once more be dropped, yielding: 



1 /A„ 
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[1 - cos(A osc s)] 



(25) 
(26) 



We see that this is just Eq. 19, but with the scale height 
H K replaced by the magnetic domain size s. For our cal- 
culations in Section III, we use the the full expressions 
(Eqs. 18 and 24) rather than the small Po approxima- 
tions. 

A final comment should be made on the topic of ALP- 
photon reconversion. Because this effect is implicit in 
the definition of Po (Eq. 11), it is accounted for in the 
photon-loss expressions Eqs. 18 and 24. However, we 
have assumed that at every height the initial ALP popu- 
lation is zero, i.e. I a {zo) = 0. Although it is not true in 
general that no axions will reach a given layer from any of 
the layers below, this is in fact a good approximation for 
our purposes. At height z , the "initial" ALP population 
is given to a good approximation by the conditions in the 
layer below, so I a (zo) = ^Po(Sz)Ij(zo — Sz). Therefore, 
Eq. 20 dictates that the reconversion rate of the initial 
axions into "new" photons is 



1 



(z) = -P (5zyi-y(z -Sz), 



(27) 



where P^rec is the intensity of reconverted photons only. 
As Po <C 1 everywhere but the resonance region, this 
quantity should be smaller than the photon-loss rate as 
long as the scale height H K (and therefore Sz) is larger 
than the width <5/i ros of the resonance region. In all cases 
of interest to us, H K ~ 100 km, whereas Sh lcs < 15 km. 



III. EFFECTS ON SOLAR ABUNDANCES 

To determine the impact of photon- ALP conversion 
in the photosphere on solar abundances, we calculated 
the expected correction factor Aloge (Eq. 7) using the 
3D hydrodynamic solar model atmosphere from [9]. We 
computed the ALP-induced opacity at wavelengths 
from 300 to 1000 nm, at each point in an (x,y,z) grid 
with resolution 50 x 50 x 500, sampled from the atmo- 
sphere simulation. Our extracted grid covers 6 x 6 Mm, 
in the horizontal directions, and optical depths in the 
range —3 < logr < 0.5, in which the majority of line for- 
mation occurs and below which conversion is negligible. 
We did this for 90 snapshots from the simulation, cor- 
responding to 45 minutes in solar time before averaging. 
All averages presented here are carried out over time and 
over the undulating surfaces of particular optical depths 
to give averages on the Tsoo-scale. We averaged the local 
values of over (x,y,t) slices of common optical depth 
at 500 nm, to arrive at expected abundance corrections 
as a function of wavelength and formation height of ar- 
bitrary spectral lines. Spatial and temporal changes in 
line profiles due to convective motions occur over the 
length scales and lifetimes of convective granules, which 
are approximately 1 Mm and lOmin for the Sun [52, 53]. 
As these scales are significantly smaller than the extent 
of the simulation, our averaged values should accurately 
represent the expected effect of ALPs on the mean solar 
spectrum. 

The continuum opacities include bound-free absorp- 
tion by the first three ionisation stages of 15 of the most 
abundant elements, as computed by the Opacity Project 
(OP) and the Iron Project (IP) [54]. The absorption by 
the anions of H, He, C, N and O is also included, as well 
as that by the molecules H 2 , H^, H^, OH, CH, H 2 
and CO - and the pressure induced absorption of tran- 
sient (Hi+Hi)-, (Hi+Hei)-, (H 2 +Hl)-, (H 2 +Hel)- and 
(H2+H2)-pairs. A complete list with detailed references 
was presented in [55]. The opacity in the 380-1644 nm 
range is dominated by the bound-free absorption by H~ , 
by free-free absorption of Hi and H - further into the 
infra-red and by bound-free absorption by Hi and vari- 
ous metals further into the ultra-violet. 

Depending on initial conditions and simulation details, 
magnetohydrodynamic models of the photosphere pre- 
dict vertical magnetic field amplitudes in the quiet Sun 
up to 0.1 T, with domains ranging from a few hundred 
to 1000 km in diameter at line-forming heights [56-58]. 
This means that Eq. 18 is valid for all our calculations; if 
s were equal to or smaller than the scale height H K , Eq. 26 
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FIG. 2. The effect of resonant photon-axion conversion on the effective line widths. Left: Abundance correction effect in 
the solar atmosphere due to particles with an axion-like coupling, as a function of optical depth r and wavelength. Here 
= 10" 4 eV and we have set g al = 0.1 GeV , to illustrate the effect, although the allowed value of g al by observational 
searches is much lower. Horizontal bars indicate the central 68% intervals of the line formation contribution functions for the 
630.0 nm [Ol] (green) and 872.7 nm [Cl] (black) forbidden neutral atomic lines. Right: horizontal cross-sections through the 
left-hand panel for (solid lines) A = 400 (magenta), 630 nm (green) and 875 nm (black). The peaks are due to resonant-photon- 
ALP conversion near optical depth log(r) = 0. Green and black dashed lines are the contribution functions dP/dlogr to the 
630.0 nm [Ol] and 872.7 nm [Cl] absorption lines, respectively. In both cases, J dP — 1 with r being the monochromatic optical 
depth at the central wavelength of the respective absorption line (rather than rsoo). 




FIG. 3. The effect of non-resonant photon-axion conversion, with a coupling large enough to solve the solar abundance problem. 
Same as Fig. 2, but with g aj = 10 GeV -1 and m a — 1 x 10 _4 eV. Solid lines on right-hand figure are for A = 400 (magenta), 
630 nm (green) and 875 nm (black); dashed lines are the Ol and Cl contribution functions (see Fig. 2 caption). The conversion 
probability remains constant throughout the line-forming region of the atmosphere, and the rise in A log e is due to the falling 
bulk opacity with height. We note that the extremely large coupling required is ruled out by several orders of magnitude. 



would be needed instead. The horizontal component of 
the magnetic field can furthermore be 2 to 5 times larger 
than the vertical component [59] . We therefore take (B) 
— 0.1 T, but caution that some studies (including [59]) 
favour a more modest 0.001 to 0.01 T. 

The impact of resonant photon-axion conversion on the 
effective line opacity is illustrated in Fig. 2 with, m a — 



4 x 10~ 5 eV and a very large value of g ai = 0.1 GeV -1 , 
giving resonant conversion at a mean optical depth of 
log(r) ~ 0.2. We also show the contribution functions for 
the 630.0 nm [Ol] and 872.7 nm [Cl] forbidden lines (re- 
produced from 3D line formation calculations presented 
in [25] and [27], respectively), indicating the contribu- 
tions from different atmospheric heights to the formation 
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of these lines. These two lines are important contributors 
to the determination of oxygen and carbon abundances 
[1-4, 9]. Due to the A^ as contribution, which may be 
seen as a suppression of the effective photon mass by the 
index of refraction of neutral hydrogen, resonant conver- 
sion within most of the line-forming region is disallowed. 
We note the sharp rise at low optical depths is because 
of the sharp drop in the opacity at these heights: Fig. 2 
b) shows the relative change in line strength, reflecting 
the decrease in continuum opacity; the absolute change 
at optical depths less than ~ — 2 has no observable effect. 

If the coupling g ai is pushed to even larger values 
<7a 7 > 10 GeV -1 , non-resonant conversion can provide 
corrections of A log e > 0.2 dex in the centre of the region 
of line-formation. We illustrate this in Figure 3. Taken 
naively, the correction factors we present here could offer 
at least a partial solution to the solar abundance problem, 
as abundances from all lines would increase relative to 
their expected values. This example combination of m a 
and g ai would induce the necessary 0.1-0.2 dex change 
in abundances from some common C and O indicators; 
similar effects would be expected with other lines used 
for abundance determination, as they typically form at 
similar optical depths. The exact amount by which the 
the abundance inferred from each line and species would 
rise would depend on line choice, line-to-line differences 
in wavelengths and formation depths, the magnetic field 
structure and so on. 

However, to this point we have not considered the via- 
bility of the couplings required, in in view of independent 
constraints from solar core or laboratory searches for ax- 
ions and ALPs. In the following Section we discuss those 
bounds and the possibility of achieving this effect with 
real particle models, showing that any observable effect 
is in fact already ruled out by existing bounds on ALPs. 

IV. DISCUSSION 
A. ALPs 

There are three types of photon- ALP conversions that 
can be probed in order to constrain g ai and m Q 3 : 

1. Conversion in hot dense cores of stars should pro- 
duce X-ray energy ALPs. Observations of horizon- 
tal branch stars constrain the rate of energy loss 
during their He-burning phase [60]. The Sun itself 
would lose significant amounts of energy from its 
core if standard ALPs with couplings large enough 
to solve the solar abundance problem truly existed, 
giving rise to observable effects in the frequencies of 
p-mode (sound wave) oscillations [61] and the solar 
neutrino flux [62]. Conversion in the Sun can also 



3 The notation M al = l/g a -y is often used in such discussions. 



be directly probed: pointing an "axion helioscope" 
at the Sun give strict limits on the reconversion of 
ALPs into X-rays [63]. The CERN Axion Solar 
Telescope (CAST) bound on Solar core ALP pro- 
duction and reconversion is g ai < 9 x 10~ n GeV -1 
for masses less than 10 -2 eV. 

2. Light shining through wall (LSW) experiments 
probe reconversion of ALPs produced from laser 
beams in an intense magnetic field. Current bounds 
are from the ALPS [64] and GammeV [65] exper- 
iments. Both shine a high-intensity 532 nm laser 
beam at a stopper in a transverse magnetic field, 
and set bounds based on the (non) observation 
of reconverted photons behind the stopper. For 
masses below 1 MeV, the ALPS bound is g ai < 
6 x 10 -8 GeV -1 whereas the GammeV bound is 
gay < 3 x 10 -7 GeV -1 . 

3. Vacuum polarisation changes from ALP conver- 
sion of a linearly polarised infra-red laser beam 
in a 2.3 T magnetic field has been probed by the 
PVLAS experiment [66]. Vacuum dichroism oc- 
curs when some of the photons oscillate into ALPs, 
causing a loss of part of the beam with polar- 
isation parallel to the magnetic field, and thus 
cause a small beam rotation. Conversely, birefrin- 
gence occurs when ALP oscillation causes a phase 
lag of one of the polarisation components, causing 
an ellipticity to develop in the beam polarisation. 
The strongest bound to date from this effect has 
been provided by the dichroism limits of PVLAS: 
9ai % 4 x 10~ 7 GeV -1 for ALP masses below ~ 1 
meV. 

These limits are much too strong to allow standard ax- 
ions or ALPs to produce an observable effect in the solar 
atmosphere. Nonetheless, there are classes of models in 
which some of these bounds do not apply. It was pointed 
out by [67] that if the ALP is a composite object anal- 
ogous to the neutral pion, then the photon-ALP conver- 
sion can be suppressed at high energies when one of the 
photons is virtual. This suppresses ALP production in 
the hot interiors of stars, alleviating constraint 1. If the 
ALPs are furthermore short-lived particles, for instance 
if they decay to a lighter particle that does not interact 
strongly with the standard model, then constraint 2 may 
also be overcome. In the latter case, the reconversion 
discussion in section III becomes irrelevant. 

However, dichroism constraints from PVLAS cannot 
be brushed aside even for non-standard models of ALPs 
that can avoid constraints 1 and 2. Indeed, the photon- 
loss that would be observed in a dichroism experiment 
is the same as the effect that would be responsible for 
photon loss in the line-forming region of the solar atmo- 
sphere. If ALPs must respect the PVLAS upper limit on 
their coupling {g ai < 10~ 7 GeV -1 ), then the coupling re- 
quired to produce an observable effect is suppressed by at 
least eight orders of magnitudes. We therefore conclude 



that current limits rule out axions and ALPs as progeni- 
tors of any observable effects on abundances derived from 
the solar photosphere. 



B. Chameleons and hidden photons 

In addition to ALPs, there exist further models of par- 
ticles which mix with the photon in an analogous manner 
to (3), and one could hope that such models may allow 
constraints to be evaded. We look at two specific models: 
chameleons, and hidden photons. 

The chameleon, first proposed as a model of dark en- 
ergy [68, 69], consists of a scalar field (j> with an effective 
coupling to matter V(4>) ~ e^' Mpl p, where /3 is a param- 
eter of the theory, Mpi is the Planck mass and p is the 
matter density. This gives rise to a particle [70-72] whose 
mass depends on the local density via oc p. Expand- 
ing the exponential coupling, the chameleon's first-order 
coupling to the photon is: 



£<t>i - -ji—F^F^i (28) 



where M^ 7 = Mpi//3 is the mass scale of new physics 
giving rise to the chameleon field. The only relevant 
phenomenological difference with ALPs will be that 
chameleons couple to photon polarisation which is or- 
thogonal, rather than parallel to, the bulk magnetic field. 
As we assume random polarisations here, this does not 
affect conversion in the solar atmosphere. The PVLAS 
experiment gives similar constraints on the chameleon 
mixing parameter as it does to the ALP case: 1/M^ 7 < 
lO^GeV- 1 ^]. 

The hidden photon (or paraphoton) paradigm [74-78] 
has its origins in an additional U{1) gauge field, with field 
strength F' which mixes very weakly with the standard 
model photon via a term 

£ 7 '7 = \xF'^. (29) 

A combination of constraints from the cosmic microwave 
background, Coulomb-interaction experiments and helio- 
scopic observations from CAST constrain the mixing pa- 
rameter x S 10~ 7 [79]. 

Crucially, the form of the coupling to photons means 
that resonant conversion of photons to either chameleons 
or hidden photons occurs in the same qualitative way as 



conversion to ALPs: 

Py^ t oc — ^ . (30) 

(m|. - < off ) + const. 

Resonant oscillation is then forbidden for either of these 
scenarios at optical depths logr < 0, due to the dom- 
inance of A^ as . Below such atmospheric heights, we 
find that the allowed parameter space of either theory 
does not permit oscillation. The hidden photon mixing 
parameter is roughly equivalent to x i — > 9ajBu/m^, 
whereas the chameleon coupling can be read directly as 
l/Mfh^, < — > g a ~f. Required parameters in both cases 
fall far above the current exclusion bounds, guarantee- 
ing that no observable conversion may take place in the 
solar atmosphere. 4 

V. CONCLUSIONS 

We have shown that for appropriate masses and cou- 
plings, particles with an axion-like interaction (Eq. 3) 
could in principle resolve the solar abundance problem, 
by inducing a slight increase in the effective contin- 
uum opacity at line-forming heights in the solar atmo- 
sphere. This in turn would reduce the computed equiv- 
alent widths of solar absorption lines for any given ele- 
mental abundance, and could bring inferred photospheric 
abundances into agreement with helioseismological re- 
sults. 

However, we found that the coupling necessary to ob- 
tain such a drastic change is ruled out by current ex- 
perimental null results. We have shown that chameleons 
and hidden photons, which exhibit a similar phenomeno- 
logical behaviour to the ALP, are also disfavoured as a 
photospheric solution to the solar abundance problem. 
We therefore conclude that new light bosons are not able 
to provide a photospheric solution to the solar abundance 
problem. 
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